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I Near-critical behavior of the free surface of an ideahy conducting hquid in an 



o 
o 



external electric field is considered. Based on an analysis of three-wave processes 
using the method of integral estimations, sufficient criteria for hard instability of a 
J>->| planar surface are formulated. It is shown that the higher-order nonlinearities do not 



saturate the instability, for which reason the growth of disturbances has an explosive 
character. 
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CIh: I. INTRODUCTION 

(N . 

^ . The electrohydrodynamic instability of the free surface of a conducting liquid in a strong 



electric field is responsible for many physical processes, such as the initiation and 

^ ' maintenance of emission of charged particles, vacuum breakdown, vacuum discharge, etc. 
^ ! The interaction of an electric field and charges induced by this field on the surface of a liquid 
' (liquid metal) leads to a growth of surface disturbances and the formation of regions with 

. a significant curvature . The dispersion law for the waves on a planar surface of an 

' ideally conducting liquid in an external electric field of strength E has the following form 



X' u^ = g\M + -\\.\'-^\\.\' (1) 

where uj is the frequency, k is the wave vector, g is the acceleration of gravity, a is the 
surface tension coefficient, and p is the density of the medium. 
It is seen from Eq. (1) that if the condition 



E"^ < Ec^ = Sny^gap, 



* Electronic address: lnick@ ami . ur an . rul 



2 



is fulfilled, then cu^ > at any |k| and, consequently, the surface disturbances do not grow 
with time. If the magnitude of the field E, which plays the role of an external governing 
parameter, exceeds a certain critical value Ec, then there arises a region of wave vectors |k|, 
for which u"^ < 0, which corresponds to an aperiodic instability. Thus, the condition E > Ec 
is a criterion for the surface instability with respect to infinitely small disturbances of the 
surface shape and of the field of velocities. 

It was shown in , where liquids with various physical properties have been consid- 

ered, that a nonlinear interaction of three standing waves that form a hexagonal structure 
can lead to a hard excitation of the instability of a charged surface. In our case, this means 
that, even at subcritical fields [E < Ec) a disturbance of a sufficient magnitude can break 
the equilibrium of a planar surface. In this connection, there arises a need of constructing 
criteria for the instability of a charged surface of a conducting liquid with respect to per- 
turbations of a finite magnitude, i.e., criteria that will permit one, proceeding from some 
initial data such as the shape of the surface and the distribution of velocities, to answer the 
question of whether or not the initial perturbation will lead to the loss of the stability of a 
planar boundary and, as a consequence, to an explosive growth of cusplike structures. This 
work is devoted to constructing such criteria using the method of integral estimations that 
was applied previously to obtain the conditions of collapse for the nonlinear Schrodinger 
equation jlol.lll|. nonlinear Klein-Gordon equation jl^-lisj]. and various modifications of the 
Boussinesq equation [l^ . 

In Section 2, we give equations of the vertex-free motion of an ideally conducting liquid 
with a free surface in an electric field and give their Hamiltonian formulation. In Section 3, 
the theory of perturbations in a small parameter, namely, in the characteristic angle of the 
surface slope, is constructed up to fourth-order terms in the Hamiltonian. The analysis of 
the surface dynamics is significantly simplified in the case of small "supercriticalities" (i.e., 
if the field only slightly exceeds the critical value Ec) 

e = {E^- El) /El, \e\ < 1, 



when only perturbations with wave numbers close to ko = \fgpJoL increase (this value of 
the wave number corresponds to the so-called dominant harmonic of surface perturbations). 
This permits us in Section 4 to construct a set of amplitude equations for describing the 
nonlinear interaction of three standing waves that form a hexagonal structure which is the 
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main interaction at near-critical values of the field E. In Section 5, we extend the method 
of integral estimations for several interacting nonlinear waves. By using this method, one 
can pass from a set of partial differential equations for complex amplitudes Ai, A2 and A3 

to to a second-order differential inequality for the norm 



by analyzing which we obtain a number of sufficient criteria of hard excitation of an electro- 
hydrodynamic instability of a charged surface. Note that most of them refer to subcritical 
values of the external electric field, when the surface is stable in a linear approximation, 
while the development of the instability is related to three- wave processes. In Section 6, we 
show, on the example of one-dimensional and square lattices of surface perturbations for 
which the three-wave interactions degenerate, that the higher-order wave processes do not 
saturate the instability but, on the contrary, lead to an explosive growth of amplitudes. 



Consider the potential motion of an ideal conducting liquid of infinite depth placed in an 
external uniform electric field of strength E. Assume that the fieldstrength vector is directed 
along the z axis and, correspondingly, in the unperturbed state the boundary of the liquid 
is a planar horizontal surface z — 0. Let the function 77(2;, y, t) specify the deviation of the 
boundary from the flatness, i.e., the region occupied by the liquid is restricted by a free 
surface z — rj. 

The velocity potential for an incompressible liquid $ satisfies Laplace's equation 




II. STARTING EQUATIONS 



= 0, 



(2) 



with the following conditions at the metal-vacuum boundary and at infinity: 



5$ (V$)^ _ {Vipf - E^ 




= -9V, 



Z^T), 



$ ^ 0, z^ -00, 
where (p is the electric-field potential. 



(3) 



4 



The first term on the right-hand side of the dynamic boundary condition (nonstationary 
Bernoulli equation) is responsible for the electrostatic pressure, the second term determines 
the capillary pressure, and the third term takes into account the effect of the gravitational 
field. The time evolution of the free surface is determined by the kinematic relationship 
(condition of non- flowing of the liquid through its boundary) 

dri 9$ „ „ ^ 
dt^~dz~ ' 

Finally, the electric-field potential ip in the absence of spatial charges satisfies Laplace's 
equation 

vV = 0, 

which should be solved simultaneously with the condition of the equipotentiality of the 
boundary of the conducting liquid and the condition of the uniformity of the field at an 
infinite distance from the surface: 

(^ = 0, z = r], 

{p —Ez, z ^ oo. 

Note that the above-written equations of motion have a Hamiltonian structure and the 
functions r]{x,y,t) and ilj{x,y,t) = ^\z=r, are canonically conjugated quantities [3] 

dip _ 6H dr] _ 6H 

'dt ~ ~Tn' dt ^ ^ 

where the Hamiltonian H coincides to an accuracy of a constant with the total energy of 
the system 

J 2 J STcp 

+/[^ + ^(v/rTTw-i)]^v. 

For a further consideration of the problem, it is convenient to represent the Hamiltonian 
in the form of a surface integral. We introduce a perturbation of the electric-field potential 
if = if + Ez. It can easily be shown that the perturbed potential if satisfies Laplace's 
equation 

VV = 0, (5) 
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with conditions 



(6) 



ip ^ 0, 2; — > 00, 



(7) 



from which it is seen that the perturbation introduced by the surface z = rj into the dis- 
tribution of the electric field decays as 2; — * 00. Taking into account that, in view of the 
incompressibility of the liquid, a relation J (p\z=rid^r = is valid, and neglecting terms whose 
variation does not contribute to the equation of motion, we obtain, using the first Green's 
formula. 



H 



I 



i/j 9$ Et] dip 
2 dn Snp dn 



ds + 



^ + ^(Vl + (Vxry)^-l 



d r. 



where ds is the surface differential and d/dn denotes the derivative in the direction of the 
normal to the surface z — rj. 

By eliminating the normal derivatives of the potentials (p and we can reduce the 
expression for the Hamiltonian to the form 



H 



I 



Er] 
Svrp 



which is more suitable for further transformations. 



III. SMALL-ANGLE APPROXIMATION 



(8) 



Our further problem is to eliminate the spatial variable z from the equations of motion, 
i.e., to pass from the initial three-dimensional equations to two-dimensional ones. To do this, 
we should write the integrand in Eq. (8) through the canonical variables 7] and ip. Then, 
there arises a need to solve Eq. (5) with conditions (6) and (7) as well as Eq. (2) with the 
condition 
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and condition (3). We use the known solutions to Laplace's equation for the half-spaces 
2; < and 2; > for functions that decay at infinity 

+ 00 +CX) 

—00 —00 

+00 +00 

1 f f z^x',y',0) 



^x,y,z) = -— / ^ ^ Tj^dx'dy', z<0. (10) 



-00 —00 



Now, we should express the magnitudes of the potentials (f and $ that enter into these 
relationships at the plane z = through their magnitudes at the boundary z — rj, i.e., 

through the functions Erj and ip. Let the characteristic angles of the surface slope be small: 
\'V±r]\ <^ 1. In this case the potentials near the z = plane can be expanded into a power 
series in surface perturbation rj: 



(^(x, y, riix, y)) = ^— -— , y, n{x, y)) = ^— —- 

z=0 n=0 



n=0 



(11) 



5=0 



By differentiating Eqs. (9) and (10) with respect to z, we find that 

where k is the two-dimensional integral operator given by the expression 

+00 +00 

kf--^[ [ ^.dx'dy'. 

27tJ J [(^x' - xy + {y' - y)^f^ 
—00 —00 

This relationship can be considered as a consequence of the fact that the Laplacian operator 
can formally be represented as 

V = (d, + (d, - k) , 

where the left-hand bracket corresponds to solutions that are asymptotically decay as 2; — > 

+00 and the right-hand bracket corresponds to solutions that decay as 2; — —00. 
By eliminating the derivative with respect to z from the expansions (11), we find 

ip{x, y, r]{x, y)) = f+(p{x, y, 0), ^{x, y, r]{x, y)) = TL$(a;, y, 0), 

where we introduced nonlinear shear operators 



n=0 n=0 
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Let ^ be the operators that are inverse with respect to the shear operators. Their form 
can be determined using the method of successive approximations 



Then, we have 



T^^ ^l^rik- + r]kr]k + ... 



These relationships, along with Eqs. (9) and (10), specify the solutions to Laplace's 
equations with necessary boundary conditions in the form of infinite series. Using these 
solutions, we can write the various possible derivatives of the potentials (p and $ that enter 
into the Hamiltonian through the functions rj and iJj. As a result, we find 



+ 



9V « 
2 p 



1 + (Vxr7)2 



To describe the initial stages of instability development on the surface of a conducting 
liquid, it is sufficient to restrict ourselves by allowance for a finite number of terms in the 
expansion of the integrands of the functional H in canonical variables. By omitting terms 
of higher than the fourth order of smallness for the surface perturbation r) and higher than 

the second order for the potential ip, (this proves to be sufficient at small supercriticalities) 
and successively integrating by parts, we finally obtain 



ipkip E'^rjkrj grf aiVviY 



(12) 



^(3) 



{Vrif - {krjf 



Tjkrjkrjkrj + rjkrf'V'^rj 



8p 



d\. 



These expressions, in combination with Eq. (4), represent a two-dimensional reduction of 
the equations of motion of a conducting liquid in an external electric field that is applicable 
if the condition of the smallness of the characteristic angles of the surface slope is fulfilled. 
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IV. AMPLITUDE EQUATIONS 

Now, let us consider the nonlinear dynamics of the perturbations of the free surface of a 
conducting hquid for the case where the magnitude of the external electric field E is close 
to its threshold value Ec, i.e., |£| <S 1. It follows from the dispersion relation (1) that, at 
small supercriticalities, only surface waves with wave numbers close to /cq can be excited. 
The main nonlinear interaction in this case will be the three-wave interaction between the 
waves whose wave vectors are turned with respect to one another by an angle of 27r/3. This 
can easily be understood from the conditions 

ki + k2 + ks = 0, |ki| = |k2| = \ks\ = ko. 

Near the threshold, it is natural to pass to envelopes using the following substitutions: 

3 
3 

where ki = {ko, 0}, k2 = {-ko/2, VSko/2}, kg = {-ko/2, -VSko/2}, and Aj and Bj 
{j — 1, 2, 3) are slow functions of the variables Xj and i/j that form orthogonal coordinate 
systems with abscissa axes directed along the wave vectors kj. Such a representation for the 

functions 1] and ijj corresponds to a hexagonal structure of the perturbed surface. 

Using these relationships for the perturbations 1] and ip we can approximate the integral 
operator k that enters into the Hamiltonian (12) by a differential operator. Let us use the 
following property: 

ke^^' = |k|e^^^ 

which is related to the fact that the Fourier transform of the operator k is equal to the 
modulus of the wave vector. Consider a plane wave of the form 

gikr _ ^i{kox+qxX+qyy) 

whose wave vector is close to ko {i.e.,\qx\ <^ ko and \qy\ <^ ko)- The quantity |k| can be 
expanded in a series in and qy-. 

|k| = y(^^Tg^)2Tg7 ^ko + qx + {2ko)-\y^ - (2A;g)-'g^?/ - (2A;o)~^g/. 
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This means that if we deal with a narrow (in the k space) wave packet with a carrying wave 
vector k = {ko, 0} (which can be represented in the form A{x, 2/)e*^°^), then the operator k 
can be approximated as follows: 



ik()X 



kA{x,y,t)e''''"' ^ {koA - iA^ - {2koy^Ayy - i(2klY^A,^yy - {2k^Y^Ayyy^ e 

(similar relations are obtained for the amplitudes Aj in the coordinates of Xj and yj). Then, 
inserting the expressions for 77 and into the Hamiltonian (12) and performing necessary 
averaging, we find (to an accuracy of terms of a higher order of smallness) 



dAj i d'^Aj 



dxj 2kQ dyj 



-3gko j (A1A2A3 + AlA*^Al)(fr. 



The dynamic equations t 
found from the relations 



lat describe the time evolution of the amplitudes A, and Bj are 
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A. B 



^* 5B*' 5A*' 
where j = 1, 2, 3. 

By varying the expression for the averaged Hamiltonian, we obtain the following equations 
for the amplitudes: 

Aj^ = koBj, 

By eliminating the amplitudes Bj from these equations and passing to dimensionless quan- 
tities using the substitutions 

r-^r/{V2ko), Aj-^A^ko, t^t/^2^o, H^2Hg/kl (13) 

we obtain the following set of equations: 

Aut = eA^ + LlA^ + 3 A^A* /2, (14) 

A2tt = eA2 + LlA2 + 'iAlAl/2, (15) 
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A^,, = eA^ + LlA^ + ?,AlAl/2, 



(16) 



where we introduced operators 



d 



I 92 



dxj dy 



2 ' 



J = 1,2,3. 



The Hamiltonian corresponding to these amphtude equations is written as follows: 

3 



H = I Y.{\AM + \L,A,\'-e\A,\ 
Li=i 



-{A^A,A, + A\AIAI) 



(17) 



Thus, we obtained equations that describe the initial stages of the development of in- 
stability of the surface of a conducting liquid in a near-critical field, when the small-angle 
approximation is valid and the main nonlinear interaction is the interaction of three standing 
waves that form a hexagonal lattice. Note that analogous equations describe the instability 
of the charged surface of liquid helium P, |l3| . 



V. CRITERION FOR EXPLOSIVE INSTABILITY 



As is known, hexagonal structures on a charged surface of various liquids are characterized 
by a hard regime of excitation 0, |^. For the set of equations (14)-(16), this means the 
possibility of an unbounded growth of the amplitude Aj in a finite time period. Indeed, in 
the simplest case, when the amplitudes are independent of the spatial variables, are real, 
and are equal to one another, i.e., Ai = A2 = A3 = A{t), the time evolution of the quantity 
A is described by an ordinary differential equation with a quadratic nonlinearity 

Att = eA + 3AV2. 

Because of its influence, the amplitude grows asymptotically {A — > 4(t — tc)^^), under 
corresponding initial conditions, i.e., the magnitude of A becomes infinite at the moment tc- 
However, what seems to be obvious for spatially uniform (coordinate-independent) solutions 
requires to be proved in the case of arbitrary amplitudes Ai , A2, and ^3. In particular, of a 
significant interest is the situation where the initial perturbation of the surface is localized 
in a certain region. 

Let us show, using the method of differential inequalities, that the nonlinear interaction of 
amplitude Aj in terms of the model (14)-(16) results in an explosive growth of perturbations 
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of the surface of a conducting liquid and find the sufficient conditions for hard excitation of 
the instabihty. To this end, we introduce the norm 



\A,\^d\ 



J = 1,2,3 



and consider the time evolution of the following nonnegative quantity: 

3 



By analogy with |12|], we doubly differentiate X with respect to t 

= J 2J2{\A,\'t + \L,A,\' + e\A/)+liA^A2As + AlAlA; 



after substituting the corresponding right-hand sides of the amplitude equations (14)-(16) 
for the multiple derivatives Aj,, and A*-^. Then, eliminating the signambiguous cubic non- 
linearity from the integrand using the expression for the Hamiltonian (17), we obtain the 
following relation: 

(18) 



Xu + 3H = -eX + [ [l^jft + l4^jf 

7 = 1 



d'r. 



Now, our problem is to approximate the right-hand side of Eq. (18) using the magnitude 
of X and thereby obtain an ordinary differential inequality. From the known Cauchy- 
Bunyakowsky integral inequality for the functions \Aj\ and \Aj\t, 



A,\]dh 



> 



it follows that J \Aj\^d'^r > Xj'^/{4:Xj). Taking also into account the obvious relations 
/ \LjAj\^d^r > for j = 1, 2, 3, we obtain from Eq. (18) 

3 2 

(19) 



Xu + 3H>-eX+^-J2^- 



Then, note that, as a consequence of the algebraic Cauchy inequality, the following rela- 
tion is valid: 



> 
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and, correspondingly, we have 

3 

Y,x,]/x,>x^/x. 

Substituting the latter inequality into (19), we obtain an ordinary differential inequality 

Xtt + 'iH>-eX + --^, (20) 

which will be the object of our consideration below. Note that analogous inequalities arise 
when deriving sufficient collapse criteria for various nonlinear partial differential equations 

The introduction of a new function Y = X~^^^ permits us to rewrite the inequality (20) 
in the form of Newton's second law 

PiY) = -l{eY^ + HY^), (21) 

where Y plays the role of the coordinate of a "particle" and P is its potential energy. 

Let the velocity of the "particle" Yt be negative (in this case Xt > 0). Then, multiplying 
(21) by Yt, we obtain 

Utit)>0, U{t)=Yt^/2 + P{Y), 

i.e., the "particle" gains an energy U upon motion. It is understandable that the sufficient 
criterion for Y to become zero and, correspondingly, for X to become infinity is the condition 
that the "particle" encounters no potential barrier even if Ut = 0, which corresponds to the 
equality sign in (2). The explosive growth of amplitudes takes place under the following 
conditions: 

(a) at £ < and if > if Y{to) < 1^1 V(3)^ and 12f/(to) < \e\l/{3H)^; 

(b) at £ < and > if 12U{to) > \e\l/{3H)^; 

(c) at £ < and F < 0; 

(d) at e > if U{to) > 0, 

In (a)-(d), t = to corresponds to the starting time moment. In this case, the moment tc, at 
which the perturbation amplitudes become infinite is estimated as follows: 

tc<to+ 



dY 







^2U{to) - 2P{Y) 
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Note that the condition l^(to) < in cases (a) and (c) is by no means necessary: after the 
reflection from a potential wall, the "particle" reaches the point Y = 0. The above conditions 
(a)-(d) can be considered as sufficient criteria of the instability of the surface of a conducting 
liquid with respect to perturbations of a finite amplitude, which distinguishes it from the 
simplest criterion of linear instability E > Ec, which was derived based on the assumption 
that the perturbations are infinitely small. Note also that conditions (a)-(c) refer to the case 
of subcritical external fields {E > Ec), when the flat surface of the conducting liquid is stable 
in the linear approximation, i.e., we deal with a hard excitation of an electrohydrodynamic 
instability. 

Thus, if conditions (a)-(d) are fulfilled, Eqs. (14)-(16) describe an infinite growth of 
amplitudes Aj. In this case, the applicability of the model (14)-(16) to the description 
of the development of an electrohydrodynamic instability is restricted by the condition of 
the smallness of the amplitudes: in the order of magnitude, the absolute values of the 
amplitudes \Aj\ should not exceed the magnitude of the parameter of supercriticality e. 
Otherwise, the model cannot be restricted to the consideration of only three-wave processes. 
As to the higher-order wave processes, there arises a question of whether they will lead 
to a stabilization of the instability or will favor an explosive growth ofperturbations (the 
experimental data of and the results of numerical calculations 0, J evidence in favor 
of the latter situation). The complexity of the estimation of their influence is related to 
the fact that the contribution of higherorder nonlinearities becomes comparable with the 
contribution of quadratic nonlinearities of the model (14)-(16), only if the amplitude of 
surface perturbations is close to the characteristic length of the wave. But in this case 
conditions of the applicability of our approach to the description of the near-critical dynamics 
of the charged surface of liquid metal based on the construction of amplitude equations 
become violated. Nevertheless, it is possible to reveal the influence of the higher-order 
nonlinearities by considering one-dimensional and square lattices of surface distributions, 
for which the three-wave interactions degenerate and the dominating interactions are the 
four-wave ones. 
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VI. FOUR- WAVE INTERACTIONS 



Let us consider perturbations of the boundary of a conducting liquid with such sym- 
metries for which the effect of three-wave processes is neghgible. Thereby, we in pure form 
separate the four- wave interactions that determine the character of the electrohydrodynamic 
instabihty at its advanced stages. 

First of all, we consider the near-critical behavior of the charged surface of a liquid metal 
in the assumption of a quasi-one-dimensional character of the arising wave. Let the wave 
vector be parallel to the abscissa axis. We pass to envelopes by using substitutions 



rj{x, y, t) = A{x, y, t)e''°^ + Ao{x, y, t)e''''^ + c.c, 

ilj{x, y, t) = B{x, y, t)e''''^ + Bo{x, y, t)e^'^'^ + c.c, 

in which the /cq ^ 2ko interaction is taken into account. Here, A, B, Aq and Bq are slowly 
varying functions of the spatial variables x and y. Substituting these relations into the 
Hamiltonian (12), we find (to an accuracy of terms of the fourth order of smallness) 



ko\B\''-2gs\A\^+-^ 



dA i d^A 



+ g\Ao\'^ + 2ko\Bo\ 



dx 2ko dy^ 
H^^^ = -2gko J [A^Ao* + A*^Ao] d\ 

The amplitude equations for the perturbations of the free surface are written in the Hamil- 
tonian form as 

dA _ 5H dB _ SH OAq _ 5H dBo _ SH 
'dt'lB^' 'dt~~5A^' ~df~SB^' ~df ~ ~5A^' 

By varying the function H, and then eliminating the quantities B and Bq, we obtain equa- 



tions of the form 



Au = 2gkoeA + 



g f d id 



a \ 2 



/co \dx 2ko dy'^ 



A + Agk^'A*A^--gko'A'A\ 



(22) 



Aqu = -2gkQAo + 4.gkQ^A^. 



(23) 
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Since the characteristic times of changes in the amphtudes at small supercriticalities are 
small {uj'^ ~ e), we neglect the derivatives with respect to time in Eq. (23). Then, the 
quantity Aq can be expressed through the amplitude A that plays the role of an order 
parameter 

Ao ^ 2koA'^. 

Using this relation, we eliminate Aq from Eq. (22) and pass to dimensionless quantities using 
scalings (13) to obtain for the complex amplitude A: 

Au = eA + LlA + sAlAl"^, s = 11/4, (24) 

to which the following expression for the Hamiltonian corresponds: 



H 



AA' + 



LiA 



2 



e\A\^ - s\A\''/2 



d\. (25) 



Note that, when neglecting the dependence of the amplitude A on y, Eq. (24) becomes 
a nonlinear Klein-Gordon equation, i.e., corresponds to the so-called model. In this 
form, it can be obtained from the equation for one-dimensional perturbations of the charged 
surface of liquid helium in the limit of the complete shielding of the field under the 
surface. Note also that if we neglect transverse modulations, then Eq. (24) coincides with 
that obtained in the Kelvin- Helmholtz theory of instability for the case of a small ratio of the 
densities of the top and bottom liquids 3]- This is due to the identity of the mathematical 
description of the planar potential flow of an incompressible liquid and a two-dimensional 
distribution of an electric field in the absence of spatial electric charges. The allowance for 
higher-order terms in the expansions in surface perturbations violates this analogy. 

Since the term \At\'^ that is responsible for the kinetic energy and the term \A\'^ that is 
responsible for the fourwave processes enter into the integrand of the Hamiltonian (25) with 
the opposite signs, Eq. (24) admits infinite solutions. This means that the cubic nonlinearity 
in Eq. (24) does not stabilize the linear instability but, on the contrary, enhances it, leading, 
under certain conditions, to an explosive growth of the amplitude A of perturbations of the 
conducting-liquid boundary. 

Another possible case when the three-wave processes are degenerate is the interaction of 
two standing waves whose wave vectors ki and \^2 are turned with respect to one another by 
an angle 7r/2 (the vector's coordinates are ki = {A;o,0} and k2 = {0, fco}). Let us represent 
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the perturbation of the surface rj in the form 

2 

+ (l2\/2 + le) ko [AiA2e'^^'+^'> + A^Ale'^^'-^'^] + c.c, 
and the perturbation of the velocity potential at the liquid boundary ip in the form 

2 

+ (6 + 4a/2) [(^1^2)* e^('^i+^2> + (A1A2*), e'C^i-i^^^j ^ ^ 

where we took into account the nonlinear interactions of the fundamental spatial harmonic 
ko with mixed harmonics 2ko and V2ko- This representation for the functions 77 and ip 
corresponds to the symmetry of a square lattice. 

Proceeding by analogy to the above-considered quasi-one-dimensional case, we obtain, 
after passing to dimcnsionless quantities, the following dynamic equations: 



Aitt 


= eAi 


+ LlAi 


+ sAi 


\Ai 


\' + <JA^ 




A2tt 


= eA2 


+ LIA2 


+ SA2 


\A2 


1^ + 0-^2 


\Ai\' 



where a = 32\/2 + 65/2, and the following designations are introduced: 

^ d i d"^ di/ i d"^ 

^ dx \/2dy'^ ^ dy ^dx^' 

The integral of motion for these equations, corresponding to the conservation of the total 
energy of a conservative system, is given by the expression 

H=J2 f {\Aj\t + \LM' - ^\Aj\' - s\Aj\y2^ dh-a I I^H^I^dV. 

The first term on the right-hand side of this functional coincides in its structure with the 
Hamiltonian (25) for the quasi-one-dimensional wave. The last term is responsible for the 
nonlinear interaction of a pair of the waves studied. Note that the coefficient cr before 
this term exceeds the coefficient s by more than an order of magnitude. This means that 
the contribution of the interaction ki <-> k2 is determining and, consequently, the square 
structure of the surface perturbations is much more favorable than the one-dimensional one. 
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In any case, for both the square and one-dimensional lattice (the latter can be considered 
as a partial case of the square lattice, corresponding to the condition A2 = 0), the four- wave 
interactions will favor the development of an instability rather than suppress it. Conditions 
for an explosive growth of the amplitudes Ai and A2 can be obtained by considering the 
evolution of the norm 

X = I + d'r. 

Acting by analogy with Section 5, we obtain the majorizing inequality 

Xu + ^H>-2eX + ^^, 



which coincides with that considered in 12|. The introduction of the variable Y = X^^/^ 
reduces the problem to the analysis of the motion of a "particle" with a coordinate F in a 
potential well P{Y) 

Analyzing this inequality for the case where the velocity of the "particle" at the initial time 
moment t = to is directed toward the origin (i.e., Yt{to) < 0), it can easily be revealed that 
the quantity Y vanishes, first, at £ > 0, if U(to) > 0, second, at £ < and H < 0, and third, 
at £ < and if > 0, if U{to) > e^/{8H) or ^^(to) < kl/(2). Here, U{Y), just as in Section 
5, denotes the total mechanical energy of the "particle." Under the above conditions, the 
norm X become infinite in a finite time, which just corresponds to an explosive growth of 
the amplitudes in the result of four-wave interactions. 

All this suggests that the higher-order nonlinearities will not suppress the explosive 
growth of amplitudes in the model (14)-(16). But in this case the above integral crite- 
ria (a)-(d) may be considered as sufficient criteria of the explosive growth of perturbations 
of the surface of a liquid metal in an external electric field. 



VII. CONCLUSION 



The main result of this work is the construction of sufficient integral criteria of instability 
for the free surface of an ideally conducting liquid in a near-critical external electric field. 
These criteria represent a generalization of the known condition of linear instability {E > Ec) 
to the case where the amplitudes of perturbations of the field of velocities and of the shape 
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of the surface are finite. The criteria found are dynamic in the sense that they take into 
account the effect of the velocity distribution in the medium at the initial time moment; 
the role of the stored kinetic energy can be decisive in the case of the hard mechanism of 
instability. 

An analysis of three-wave and four-wave nonlinear interactions (this corresponds to the 
allowance for quadratic and cubic nonlinearities in the amplitude equations) showed that the 
development of the electrohydrodynamic instability has an explosive character, i.e., leads to 
the appearance of singularities in the solutions in a finite time. This conclusion qualitatively 
agrees with the results of numerical simulation of the development of the instability of the 
boundary of liquid metal: it was shown in jj] that the curvature of the surface increases 
according to a power law characteristic of the explosive instability and causes the formation 
of specific features of a cusplike type. 

Note in conclusion that the criteria of hard instability analogous to the criteria (a)-(d) 
can also be obtained for insulating liquids with induced surface charges j^, for insulating 



liquids with free sur 



to this category) 



ace charges (liquid helium and liquid hydrogen in an electric field refer 



and for ferromagnetic liquids in a vertical magnetic field. 
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